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Abstract 
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1 Introduction 



In the path integral formalism, it is well-known that the semiclassical or the WKB approximation 
is a straightforward and moreover the most plausible prescription in nonpertubative treatments 0. 
Under that, first we must find out solution(s) of the equation deriving from the condition that the 
first (functional) derivative of an action, a quantity in the exponent of path integration, should 
vanish. Second we should check a positivity of the second derivative of the action at that point (or 
trajectory), that is, convexity of the action. If a solution meets this condition, designated as the 
stability condition it is called the saddle point. Third, expanding the action around the saddle 
point and regarding terms greater than the third order derivative as perturbation we have a series 
expansion expressed by a loop expansion parameter, such as ft or 1/N. Almost every cases follow 
these standard procedures but sometimes there happens a situation that the second derivative itself 
also vanishes at the point, which is so called a caustics. The cure for this case has already been given, 
for example, in the textbook by SchulmanJ^: we should start with a point given on the condition 
that the second derivative vanishes. We then utilize the Airy integration, instead of the Gaussian 
integration in the standard case. A thorough study for caustics (quantum caustics) is seen in the 
recent works by Horie et al.|lj. 

Another powerful technique exclusively used in the path integral formalism is the auxiliary field 
method(AFM); characterized by an additional degree inserted into the original expression in terms of 
a Gaussian integral to kill a four-body interaction. Although the method has a long history|S]0, it has 
recently been clarified that we can improve results when going into higher orders of the loop expansion 
in terms of the auxiliary field. It is important to note that even if the loop expansion parameter is 
not tiny but unity, results are still persistent]?]. There are many approximation schemes even under 
the path integral formalism, such as variational method|Sj, optimized perturbation theory^, dilute 
gas approximation lOj, valley method [TT]. etc., but "Simpler is Better". 

In this paper, we study a 4-fermi interaction in O-dimension as a prototype of the Nambu-Jona- 
Lasinio(NJL) model, 12 . Although the NJL model was originally handled by the mean field technique, 
it is extensively analyzed with the aid of AFMpJ- Our motivation is to establish the validity of AFM 
in NJL as was in the bosonic case|7j, however, we have found that the model itself posseses caustics. 
Therefore it is tempting to analyze it as a toy model to check the powerfulness of AFM even in a 
caustic example. By employing the standard prescription for the caustics, it is understood that AFM 
can give us a more accurate result if higher orders are taken into account. 

The paper is organized as follows: in §2, we propose our model and calculate its integration. We 
then introduce an auxiliary field and apply the saddle point method to that. In §3 and §4, we study 
cases of the 4-fermi coupling, A 2 , being posive and negative, respectively. The negative case, in §4, 
contains caustics then we follow the standard procedure, utilizing the Airy integration instead of the 
Gaussian one. The final section is devoted to discussions. 



2 Grassmann Integral and Auxiliary Field Method 
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The model considered here is 
Z = 
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and the coupling constant A 2 is supposed real. We have introduced 2iV-Grassmann variables and our 
notaion is followed from the textbook of ref.[T). 

Z in eq.Q is exactly calculable by expanding the exponent and using the binomial theorem: 
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Recalling the formula, 



we obtain 
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with [q] being the Gauss' symbol. In the following we compare this result with that of AFM with 
the help of the loop expansion, whose expansion parameter is now 1/N. 

Introduce an auxiliary field, y, to cancel the (£* • £) 2 term in eq.Q: insert the identity 
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with 
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When AT goes larger (although we will assume that JV is not so large in the following), it is 
expected that the integral is dominated by the saddle point yo obeying both conditions such that 



= S m (y ) = y - 
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and 



S (Vo) > ; stability condition 



Here and hereafter we employ the notation, 

S (n) (y) 

Write 



d n S(y) 



dy" 

O = lo + Xy , 



to find that eq.JSJ becomes 
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Expanding S(y) around yo, 
ay exp 



then making a change of variable, y — yo )—* y/y/N, we obtain 
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with = S (n) (2/o)- In view of eg. (1141 the terms of n > 3 becomes ineffective because of the factor 
w j 1 j c j 1 ena bi e us to regard these as perturbation. (As stressed before, however, N will be set 
in a small finite value in the following.) The first 6 terms are expressed as 
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respectively. The expansion up to the 3-loop correction is therefore read as 
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3 Calculation: A 2 > 

In this case, ea. 11311 reads 

ft± 



1± VT + 4A 2 



so that the stability condition, eq.@, with eg. lib! gives 
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that is, both solutions meet the stability condition thus should be taken as the saddle points. (We 
have put u> — 1 here and hereafter.) Therefore from ea. l|21jl the L-loop approximation for Z is given 

by, 
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Figure 1: The ratio; In Z tTee / lnZ exact for 1 < A < 10. The solid, the dotted and the dashed-dotted lines 
desiganate N — 2, 10, and 100 cases, respectively. When N becomes larger the ratio closes to unity, showing 
reliablity of large N approximation. 
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with the aid of eqs. H5ll ~ 120l l. In eq. 12611 we have the 1- or the 2-loop approximation if we neglect 
0(1/N) or C(l/iV 2 ) term respectively. 

To see how the results depend on N, we plot the ratio, lnZt roo to lnZ oxac t in Fig^for the cases 
of iV = 2, 10 and 100 and 1 < A < 10. (We should need the logarithm since Z itself becomes huge in 
N = 100.) 

From Fig^ it is obvious that the tree result becomes exact under iV i— > oo, which guarantees the 
large iV expansion. However, as is the bosonic case|7j, even when TV is small, the smallest, N — 2, 
in this case, we can get a better answer by taking higher loops into account: the results are listed in 
Table where the ratio ZL-ioop/Zcxact is listed for L — 0, 1, 2, 3. In spite of the fact that at A = 10 
the error is reaching ~ 50% in the tree approximation, it is cured to 4% in the 1-loop correction and 
almost to the exact value in the 2-loop. The same ratio is plotted in Figs|5| which shows that higher 
the loop we go the better the results are. As was discussed in the previous work|7j, however, since 
the loop expansion is merely an asymptotic expansion, it is not always true that a higher loop effect 
can be more accurate; whose fact can be seen in Fig[jJ] The 2-loop result becomes better than that 
of the 3-loop in 2 < A < 4. 
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A 


exact 


tree 
tree/ ex. 


1-loop 
1-loop/ ex. 


2-loop 
2-loop / ex. 


3-loop 
3-loop / ex. 


icr 3 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.00 


1.00 


1.00 


1.00 


10" 2 


1.0001 


1.0001 


1.0001 


1.0001 


1.0001 


1.00 


1.00 


1.00 


1.00 


10- 1 


1.0050 


1.0100 


1.0050 


1.0050 


1.0050 


1.00 


1.00 


1.00 


1.00 


1 


1.5000 


1.8147 


1.5346 


1.5039 


1.4996 


1.21 


1.02 


1.00 


1.00 


10 


51.000 


74.679 


53.050 


50.866 


50.902 


1.46 


1.04 


0.997 


0.998 



Table 1: The exact value as well as that of Zl-i o P and the ratio, ^L-ioop/^cxact, are listed in N = 2 case for 
10 -3 < A < 10. The 2-loop results almost reproduce the exact value. 
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1- loop/exact 
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3- loop/exact 
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1- loop/exact — 

2- loop/exact 

3- loop/exact 



Figure 2: The ratio; Z^^i oop /Z cxac t in N = 2 case for 10 3 < A < 10. The solid, dotted, dashed-dotted and 
dashed-double-dotted lines correspond to the tree, 1-, 2-, and 3-loop results, respectively. 



4 Calculation: A 2 < ; caustic case 

Next, we consider the A 2 < case. Write A 2 = -|A| 2 to find Q and S^ 2) are 



A± = V o — . (27) 



respectively. From ea. i'IHl . we find Sq 2 ' = when |A| = 1/2, that is, caustics. We cannot proceed 
anymore around this value with the standard recipe in terms of Gaussian integration. We then divide 
the region of |A| into three segments: i) |A| < 1/2 , ii) |A| > 1/2, and iii) |A| ~ 1/2, and designate 
i) and ii) as the Gaussian region and iii) as the Airy region. The reason is that in i) and ii) we can 
proceed as usual but in iii) we should adopt the Airy integral, as will be shown in the following. 
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Figure 3: The 2- and 3-loop results to the exact value are plotted. The 2-loop dominates the 3-loop between 
A rsj 2 and A ~ 4, implying that the loop expansion is merely an aymptotic expansion. 

4.1 Gaussian region I: |A| < 1/2 

Let us start with the case, |A| < 1/2. In view of eg . 112811 . only i?+ is the saddle point; since Q- does 
not satisfy the stability condition. A general expression, eg. 1 1211 . then reads 
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4.2 Gaussian region II: |A| > 1/2 

In this case, eas. 12711 . 128H are rewritten as 
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As Sg is complex, the integration range should be analytically continued from the real axis to the 
contour satisfied with conditions, 
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with 



V4|A| 2 - 1 
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and 
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The one of the conditions 13211 thus reads 
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and the rest of it gives 
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which is apparently positive. Therefore the stability condition is established for both cases i?+ and 
S7-. The quantities, eos.<l3TJl.ltT51. and CT ~ read 
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As was in eg. (1241 both solutions Q± come into play but they are complex conjugate each other, so 
that ea. 1241 1 reads 



^L-loop = Z L-loop + Z L-loop = 2Re 
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(48) 



where use has been made of the Gaussian integration formula, 
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Here in eg. 1481 . Z^oop anc ^ ^2-loop approximations are obtained by dropping terms of 0(1/N) and 
0(1/N 2 ) respectively. 

4.3 Airy region: |A| 1/2 

As was mentioned several times before, the Gaussian integtration cannot be used in this region 
because w 0; instead, we should employ a new approximation based on the Airy integral: first 
find the solution y of 
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We should adopt Q + as a suitable solution; because which connects smoothly with the saddle point: 
D+ yields 5^=0. At this point 
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Second, expand the exponent S(y) around y and make a change of variable, 
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to find 
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where 
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Note a ~ 0(1) even under AT i— > oo since f — 2|A| « in this region. Therefore, third we keep the 
first two terms in the exponent and regard others as perturbation to obtain 
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Figure 4: iV = 2 case for 10 -3 < |A| < 1: the solid, dotted, and dashed-dotted lines correspond to the exact, 
the Gaussian 3-loop, the (Airy) 3-loop results, respectively. The Airy curve drops down at |A| ~ 3.5, reaches the 
bottom ~ —1.88 at |A| = 0.238 and goes up to infinity as |A| i— ► 0. 



This expansion is not a loop expansion given in terms of 1/N but of 1/N 1/I3 , however, we adopt 
the same terminology of the loop-expansion as the Gaussian case, that is, designate the first factor 
as a tree, the second integral of 0(1) as a 1-loop, and those of 0(1/ 'N 1 * 13 ) and 0(l/N 2 ^ 3 ) as 2- and 
3-loop, respectively. We thus write ^3_i 00 p as an approximation keeping terms up to 0(l/N 2 ^ 3 ) in 
ea. l|61|l . Finally, we utilize the Airy integral [13] 
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Im[a] = 
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Figure 5: The same plot for N — 10. The 3-loop Airy at 0.4 < |A| < 0.6 as well as the 3-loop Gaussian at 
|A| < 0.4 and |A| > 0.6 fits excellently everywhere, showing reliablity of large N approximation. 



|A| 


exact 


tree 
tree/ ex. 


1-loop 
1-loop/ ex. 


2-loop 
2-loop / ex. 


3-loop 
3-loop / ex. 


0.3 


0.955 


20.832 


1.7111 


1.2467 


0.7789 


21.8 


1.79 


1.31 


0.816 


0.4 


0.92 


1.5180 


0.7305 


0.9566 


0.9341 


1.65 


0.794 


1.04 


1.02 


0.5 


0.875 


0.6796 


0.6788 


0.8752 


0.8752 


0.777 


0.776 


1.00 


1.00 


0.6 


0.82 


0.5615 


0.7608 


0.8419 


0.8145 


0.684 


0.928 


1.03 


0.993 


0.7 


0.755 


0.5888 


0.8820 


0.7758 


0.7411 


0.780 


1.17 


1.03 


0.982 


0.8 


0.68 


0.6813 


1.0105 


0.6494 


0.6628 


1.00 


1.49 


0.955 


0.975 



Table 2: The Airy loop expansion in N = 2 case for 0.3 < |A| < 0.8. As is the case of the Gaussian, the 3-loop 
correction improves the results but in 0.4 < |A| < 0.7 the 2-loop correction is sufficient. 



with Ai(a) being the Airy function, instead of the Gaussian integral. Therefore 
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(64) 



In ea. (l63H . if we drop 0(1) (Ai(a)) term, we have the tree result, and if drop O^l/N 1 ^ 3 ) term, we 
have the 1-loop result and so on. By noting that 



Ai (2) (a) = aAi(a) , 



(65) 



which is immediately derived from the definition, eq. (l62l l. the expression, eg. 1631 . can be given by 
the Airy function and its first derivative. After that we estimate it numerically using Mathematica. 

In order to show the TV-dependence of the results, we first plot the whole region, 10" 3 < |A| < 1, 
in Fig|l]for N — 2 and in Fig|S|for N = 10. As is seen from the graphs, fitting of TV = 10 in Fig|5| 
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|A| 


exact 


Gauss 
Gauss/ ex. 


Airy 
Airy/ ex. 


|A| 


exact 


Gauss 
Gauss/ ex. 


Airy 
Airy/ ex. 


10 -3 


1.0000 


1.0000 
1.00 


— 


0.5 


0.8750 


CO 


0.8752 
1.00 


10~ 2 


1.0000 


1.0000 
1.00 


— 


0.6 


0.8200 


0.3577 
0.436 


0.8145 
0.993 


ict 1 


0.9950 


0.9950 
1.00 




0.7 


0.7550 


0.6764 
0.896 


0.7411 
0.982 


0.2 


0.9800 


0.9800 
1.00 




0.8 


0.6800 


0.6500 
0.956 


0.6628 
0.975 


0.3 


0.9550 


0.9556 
1.00 


0.7789 
0.816 


0.9 


0.5950 


0.5793 
0.974 


0.5941 
0.999 


0.4 


0.9200 


0.9614 
1.04 


0.9341 
1.02 


1.0 


0.5000 


0.4904 
0.981 


0.5486 
1.10 



Table 3: Results in N = 2: those obtained by the Gaussian 3-loop expansion are listed in 10 -3 < |A| < 1 except 
|A| 7^ 1/2. Also those by the Airy 3-loop in 0.3 < |A| < 1. The lower quantities correspond to the ratio of Gaussian 
and Airy results to the exact value. 

becomes better than that of N = 2 in Fig0] so that in the following we deal only with N — 2 to 
check the reliablity of AFM. It is also apparent that the Gaussian loop expansion matches with the 
exact value for |A| > 1 thus here and hereafter we concentrate only on the region, 10" 3 < |Aj < l 1 . 

To see efficiency of the Airy loop expansion as in the Gaussian case, A 2 > , of Table □ we 
list the Airy L-loop approximation for 0.3 < |A| < 0.8 in Table |2] From this, although the 3-loop 
approximation does work well for 0.4 < |A| < 0.8, the Airy 2-loop correction is enough to produce an 
accurate value for 0.4 < |A| < 0.7. The situation is similar in A 2 > 0. Therefore the final result with 
the Airy 3-loop in N = 2 is seen in Table |3] From this, if we incoperate the Airy region by switching 
back and forth, from and to the Gaussian region, around jAj ~ 0.4 and ~ 0.7, we find that errors 
remain within 2%. 

5 Discussion 

In this study, we show that AFM for a zero-dimensional 4-fermi model does work well as in the bosonic 
case|7| when taking higher "loop" — 1/7V expansion in the Gaussian case and 1/N 1 ' 3 expansion in 
the Airy case — effects into account even if N is in the smallest number N = 2. 

So far we have assumed that A 2 is real but of course it is not necessary so. However, there arises 
nothing new when we generalize A 2 to be complex. In the complex A plane, there are two areas, located 
on the imaginary axis around A ~ ±i/2 (see eg. 12311 i. where usual saddle point approximation breaks 
down; caustics. Other regions can be treated in a usual manner, that is, more accurate answer can 
be obtained when taking higher- loop effects into account. 

Armed with this study, we should proceed to check the validity of AFM in quantum mechanical 
models. The work is now under progress|14|. After this, if we would ascertain the validity of AFM 
further, we could draw a conclusion of the issue; whether the quantum effect of pions does wipe out 
the symmetry breaking 15] or not|16|. 
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